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SUMMARY 


Preflight thermal stress analysis of the space shuttle orbiter wing skin panel 
and the thermal protection system (TPS) was performed. The heated skin panel ana- 
lyzed was rectangular in shape and contained a small square cool region at its 
center. The wing skin immediately outside the cool region was found to be close 
to the state of elastic instability in the chordwise direction based on the con- 
servative temperature distribution. The wing skin was found to be quite stable 
in the spanwise direction. The potential wing skin thermal instability was not 
severe enough to tear apart the strain isolation pad (SIP) layer. 

Also , the preflight thermal stress analysis was performed on the TPS tile under 
the most severe temperature gradient during the simulated reentry heating. The 
tensile thermal stress induced in the TPS tile was found to be much lower than the 
tensile strength of the TPS material. The thermal bending of the TPS tile was not 
severe enough to cause tearing of the SIP layer. 


INTRODUCTION 


The space shuttle orbiter reenters the earth's atmosphere at extremely high 
velocity (almost Mach 25 at the reentry altitude of 400,000 ft). Because of severe 
reentry aerodynamic heating encountered by the shuttle orbiter, the entire shuttle 
structure is covered with a thermal protection system (TPS). 

Three types of TPS are used to protect the orbiter: (1) felt reusable surface 

insulation (FRSI), (2) low- temperature reusable surface insulation (LRSI), and (3) 
hiqh-temperature reusable surface insulation (HRSI). The TRSI is a highly flexible 
blanket-type layer of heat shield, which is bonded to the low-heating zones of the 
orbiter (for example, the wing upper surface and the outer surface of the payload 
bay door) with room-temperature vulcanized (RTV ) rubber. The LRSI (white tiles) and 
HRSI (black tiles) are the low-density porous silica tiles. Each tile, which has an 
average size of 6 in by 6 in, is individually bonded with RTV rubber to the orbiter 
structure with a strain isolation pad (SIP) layer sandwiched in between. The SIP 
layer is highly flexible and serves as a cushion to absorb the thermal strain 
incompatibility between the TPS tiles and the aluminum skins. The LRSI is used 
in the medium heating zones, such as the upper region of the glove and the wing 
upper surface near the leading edges. The HRSI is used in the highly heated areas, 
such as the lower surfaces of the wings and fuselage. Some of the gaps between the 
TPS tiles in the highly heated areas are filled with ceramic-coated alumina mats 
(gap fillers) to prevent hot gases from coming in contact with the substructure at 
the bottom of the gaps. 

Before the first space transportation system (STS) flight of the space shuttle 
orbiter Columbia, all possible critical structural problems induced by the reentry 
heating had to be investigated. One concern was the possible thermal buckling of 
the orbiter skin panels. Although the SIP layer is designed to absorb the possible 
skin buckling effect on the TPS tiles, in case the skin buckling is too severe, 
the TPS tiles miqht debond from the skin, causing a partial or total loss of the 
TPS function. 


Another concern before the STS flights was the possible loss of the TPS tiles 
during the forced convective cooling period after the end of the reentry heating 
when the structural temperatures have reached their peak values* At this time, if 
some TPS tiles should detach from the skin and cause the skin to be exposed to the 
outside cool air, the thermal gradient buildup in the skin might cause an elastic 
instability in this region and cause more TPS tiles to detach from the skin. 

Another key concern was the effect of TPS thermal deformation on the SIP layer. 
During the reentry when the TPS is under the most severe thermal gradient, the 
resulting TPS thermal deformation might pull the SIP layer apart and cause debond- 
ing of the TPS. 

The purpose of this report is to perform the analysis on these potential problems 
and to draw some preliminary conclusions on the severity of their effects. 
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DESCRIPTION OF WING SKIN PANEL 


Figure 1 shows the location of the space shuttle orbiter left wing lower skin 
panel that was investigated. The rectangular skin panel is bounded by stations 
x 0 1 249 and x G 1 307 in the chordwise direction, and by stations y 0 -226 and y 0 -254 

in the spanwise direction. The reinforcing hat stringers are oriented ii} the 
spanwise direction. Figure 2 shows an enlarged view of the rectangular panel. 
Consider that the panel has been heated up to a peak temperature of 320°F (the 
estimated worst case of heating) near the end of reentry. Suppose one TPS tile 
(6 in by 6 in) located in the center region of the panel is lost after the simu- 
lated reentry heating and causes the wing skin in the central region of the 6-in 
by 6-in area to be cooled down to some lower temperature, say, 180°F, by the out- 
side cool air. The resulting thermal gradient buildup might induce thermal buck- 
ling in the periphery of the cooled region (fig. 2). The problem is to calculate 
the induced thermal stresses in the panel, to estimate the degree of severity of 
the possible thermal buckling, and to explore if the skin buckling will cause the 
TPS tiles adjacent to the cooled region to debond from the wing skin. 


THERMAL STRESS ANALYSIS 


A NASA structural analysis (NASTRAN) finite-element computer program (refs. 1 
and 2) was used in the thermal stress analysis of the orbiter wing skin panel. In 
setting up the NASTRAN structural model (fig. 3), the hat-stringer reinforced skin 
panel was replaced with an equivalent panel of uniform effective thickness. The 
model has 117 grid points and 96 quadrilateral membrane and bending elements 
(CQUAD2, refs. 1 and 2). The cool and hot regions were maintained at constant 
temperatures of 180°F and 320°F, respectively. The boundary conditions used were 
(1) free edges and (2) fixed edges. 

Figure 4 shows distributions of thermal stress Gy(x) along the x-axis predicted 

from the NASTRAN. The solid stepped curve is for the free edges; the dashed stepped 
curve is for the fixed edges. The peak compressive stress near the boundary of the 
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cool region for the case of free edges is o y (x) = 5140 lb/in 2 . The cool region is 

under considerable tension. If the NASTRAN elements were refined, the value of the 
peak compressive stress might become larger. Figure 5 shows the distributions of 
thermal stress o x (y) along the y-axis calculated from the NASTRAN. The stepped 
solid and dashed curves are for the free and the fixed edges, respectively. For 
the free-edge condition, the NASTRAN gives o x (y) = 4960 lb/in 2 for the peak com- 
pressive stress near the boundary of the cool region. Again, refinement of the 
NASTRAN elements may increase the value of the peak compressive stress. 


BUCKLING ANALYSIS 


After the peak values of the compressive stresses near the boundary of the cool 
region are known, the next problem to consider is the possible buckling of the skin 
panel element immediately outside the cool region. 

Consider a small rectangular skin element 1 (fig. 6) of length p and width i, 
located between the two adjacent hat stringers and in the highly compressed region 
(the immediate exterior of the cool region). For simplicity, the deformed shape of 
the rectangular element may be assumed to be 


h . M A ■ » Y / - \ 

sm — sxn -ir- ( 1 ) 

P * 

where w is the deflection of the rectangular plate element and h is its peak value 
at the center of the rectangular plate element. Based on this deformed shape, the 
buckling load (N x ) cr for a simply supported rectangular plate uniformly compressed 

in the x-direction may be expressed as (taking m = n = 1, ref. 3) 


(N x ) 


x'cr 


? (- • SI 


( 2 ) 


where m and n are the numbers of half-waves in the x- and y-directions , respec- 
tively, and D is given by 


E*t3 


D = 


12(1 - vt ) 


(3) 


where E a is the modulus of elasticity of the aluminum material and t is the thick- 
ness of the element. Combining equations (2) and (3), the buckling stress (o x ) cr 
may be expressed as 


(°x)cr 


t 


ir 2 E 


a 




where v a is Poisson's ratio of the aluminum material. 


(4) 
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For the shuttle skin, t = 0,062 in and p = 1.75 in. Taking E a = 10 x 1 0 6 lb/in 2 , 
v a = 0.33, and Z - 2 in, equation (4) yields 

(o x)cr * 3612 lb/in 2 (5) 

This value lies between the two peak compressive stresses o x = 2650 lb/in 2 for the 
fixed-edge condition and a x = 4960 lb/in 2 for the free-edge condition (fig. 5). If 

the actual edge condition is close to the free-edge condition, then the skin panel 
immediately outside the cool region is close to an elastic instability in the 
x-direction . 

Next, consider the second rectangular skin element 2 (fig. 6), located between 
the two adjacent hat stringers next to the boundary of the cool region. For this 
case, the length of the plate is 6 in and the width is 1.75 in. Using equation (4) 
and replacing (o x ) cr , p, and Z, respectively, with (o y ) cr , 6 in, and 1.75 in, then 
(0y) cr *n a Y he calculated as 


(a y ) cr = 169,340 lb/in 2 (6) 

This value is many times greater than the peak compressive stress a y predicted from 

the NASTRAN (fig. 4) and indicates that no buckling will occur in the direction of 
the hat stringers. 


Deflection Curve 

As shown above, the wing skin near the boundary of the cool region could buckle 
in the x-direction (chordwise direction). The next problem is to calculate the 
panel peak deflection h (fig. 6) and to determine if h will cause debonding between 
the TPS tile and the SIP layer. 

Based on equation (1), the deflection curve of the rectangular skin element in 
£ 

the y = -j plane may be written as (fig. 6) 

w = h sin — (7 ) 

P 

From equation (7), the slope of the deflection curve is given by 


dw 

dx 


irh 


cos 
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The small element of the curve length dp along the above deflection curve may 
be expressed as 


dp = 



(9) 


6 


Combining equations (8) and (9), the total curve length p of the deflection curve 
may be expressed as 


P = 



dx 


which may be expressed in terms of an elliptic function as (see appendix) 


P 
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t F FF -(*• f) 


( 10 ) 


( 11 ) 


where E(k, -?j-) i s complete elliptic integral of the second kind, and k is the 

modulus of the elliptic function given by 


k 



( 12 ) 


As h << p, the integrand of equation (10) may be expanded and integrated as 
(see appendix) 


I 


I 



from which h may be expressed as 



(13) 


(14) 


(15) 


Taking p = 1.75 in, E a = 10 x 10 6 lb/in 2 , and o x * 4960 lb/in 2 (the peak com- 
pressive stresses predicted from NASTRAN for free edges, shown in figure 3), the 
value of h may be calculated from equation (15) to yield 


h = 0.0248 in 


(16) 


Stress Induced in SIP 

After the deformation of the skin panel is known, the next problem to consider 
is whether the above skin deformation h could overstress the SIP layer. The thick- 
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ness t s of the SIP layer is 0.16 in. Thus, the maximum strain e s in the SIP thick- 
ness direction induced by the skin buckling is 


h 0.0248 

S tg 0.16 


0.1550 in/in 


(17) 


From reference 4, the modulus of elasticity ]E S for the SIP material is E s * 

117 lb/in 2 . Therefore, the induced stress a s in the SIP thickness direction is 


°s = E s £g = 117 X 0.1550 = 18.14 lb/in 2 


(18) 


This value is considerably lower than 40 lb/in 2 , the tensile strength of the SIP 
material (ref. 4) s It is therefore unlikely that the skin buckling will cause 
tearing of the SIP layer. The next problem that must be considered is the thermal 
interaction between the TPS and the SIP layer. 


THERMAL INTERACTION BETWEEN TPS AND SIP 


When the outer surface of the TPS is heated to high temperature during reentry, 
the TPS will deform due to a severe thermal gradient across its thickness. The 
problem is to calculate the thermal stresses generated in the TPS and to determine 
whether or not the TPS deformation will tear apart the SIP layer. 


Thermal Stresses in TPS 

Figure 7 shows the temperature distribution across the thickness of a TPS tile 
at 450 sec after simulated reentry, when the temperature gradient reaches the 
maximum value. The temperature curve shown in figure 7 was obtained from a one- 
dimensional heat-transfer analysis using simulated aerodynamic heating for the STS-1 
flight trajectory (refs. 5 to 8). This temperature distribution will cause the TPS 
to deform and induce a stress distribution in the thickness direction of the SIP, 
as shown in figure 8. Assuming that the temperature distribution in the TPS is 
a function of z only (fig. 7) and neglecting the low shear resistance of SIP, the 
thermal stresses in the TPS in the x- and y-directions, o x and a y , may be written as 

b I^T 

a x (z) = a y (z) = - —r—— [t(z ) - T a ] (19) 

where i s the coefficient of thermal expansion of the TPS material, E^ is the * 

modulus of elasticity of the TPS, v T is Poisson's ratio of the TPS, the quantity 

a T E T I 

is associated with biaxial deformation of the TPS tile in the x- and y- 

1 - v T 
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directions, T(z) is the temperature at any point z in the TPS tile, and T a is the 
average TPS temperature across its thickness (fig. 8 ) and is given by 


-w: 


T(z) dz 


( 20 ) 


where 2c is the thickness of the TPS tile. 

The temperature curve shown in figure 7 may be represented by the following 
quadratic function of z: 

T(z) = a^z 2 + a 2 z + a 3 (21) 

where the coefficients a-| , a 2 , and a 3 are given by ai = 609. 58°F/in 2 , a 2 = 946°F/in, 
and a 3 = 446.42°F. 

Based on the temperature profile given in equation (21), equation (20) may be 
integrated as 


Substitution of equations 

a x (z) 


a.c - 3 

T a = - 3 — + a 3 

(21) and (22) into equation (19) yields 


= o v (z) = - 


Vt 

1 " Vn 


&i K z2 ■ + a2z ] 


( 22 ) 


(23) 


The maximum tensile stress is located at the lower surface of the TPS tile and is 
given by setting z = -c in equation (23); hence, 

Vt ( 2 a i c2 \ 

o x ’(-c) = Oy(-c) = 1 _ v - \a 2 c - — ~ — j (24) 

For the TPS with a 9-lb/ft3 weight density, the physical properties and dimen- 
sions (ref. 7) are as follows: the modulus of elasticity Et = 2.5 x 10 ^ lb/in2, 

-7 

Poisson’s ratio = 0.16, the coefficient of thermal expansion a T = 3.3684 x 10 

in/in-°F, the tensile strength Of = 21 lb/in 2 * the thickness 2c = 2.54 in, and the 

length 2L = 6 in. Using these data and the values given previously for a-| and a 2 / 
equation (24) gives the maximum tensile stress in the TPS as 

0 x (-c) = Oy(-c) = 4.02 lb/in 2 (25) 

This value is only 19 percent of the TPS tensile strength Of = 21 lb/in 2 . Thus, 
the TPS thermal deformation will not induce tensile failure at its lower surface. 


Stresses Induced in SIP 


The bending moment M T induced by the thermal stress a x (z) in the TPS is given by 


Mrp = 



0 x (z)z dz 


(26) 


Substituting equation (23) into equation (24) and integrating results in 


Mrp = — 



(27) 


This moment M T must be balanced by the resisting moment generated by the nonuniform 
distribution of the normal stress 0 2 in the SIP thickness direction (fig. 8). The 

center region of the SIP will be under tension , and the regions near the edge of the 
SIP will be under compression because the TPS will bow upward under reentry heating. 

Let o z be expressed as a quadratic function of x: 

0 z (x) = b<| x 2 + b 2 x + b 3 (28) 

From figure 8, the coefficients r b 2 , and b 3 can be determined as 


= - 

°t + o c 

L2 

(29) 

b 2 

= 0 

(30) 

b 3 

* a t 

(31 ) 


where 0 t and o c are, respectively, the magnitudes of the peak tensile stress and 

the peak compressive stress in the SIP thickness direction. Substituting equa- 
tions (29) to (31) into equation (28) yields 


x 2 

0 2 (x) -* - (o t + °c^rj + °t 


(32) 


The two stresses a t and o c may now be determined from the balance of forces and 
the balance of moments. For the balance of forces. 


For the balance of moments. 



0 2 (x) dx = 0 


0 z (x)x dx = Mt 


(33) 


(34) 
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From equations (32) to (34), a t and o c may be expressed as 



°c - 2 0 t 

Substituting equation (27) into equation (35) yields 


8 


0 4 - = ~ 


Ot Em 

IJI x 


3 \ 1 - V„ 


a2 c ~ 

L 2 


Using the previously given known data, equation (37) yields 


(35) 

(36) 


(37) 


o t = 15.58 lb/in 2 (38) 

t 

it the aluminum skin deforms in the opposite direction of the TPS (worst case), 
then the total tensile stress induced in the SIP will be 


o s + o t = 18.14 + 15.58 = 33.72 lb/in 2 (39) 

This value is 85 percent of the SIP tensile strength, 40 lb/in 2 (ref. 7). Thus, it 
is unlikely that the combined deformations of the aluminum skin panel and TPS will 
tear the SIP layer. 


CONCLUSIONS 


Preflight thermal stress analysis was performed on a heated space shuttle wing 
rectangular skin panel containing a square cool region at its center. The analysis 
showed that based on a conservative temperature distribution, the wing skin imme- 
diately outside the cool region was close to an elastic instability in the direc- 
tion transverse to the hat stringers. However, the wing skin thermal deformation 
was found not to be severe enough to cause tearing of the strain isolation pad 
(SIP) layer. 

A preflight thermal stress analysis was also conducted on the thermal protec- 
tion system (TPS) tile subjected to the most severe simulated reentry heating. 

Based on this analysis, the thermal stress generated in the TPS tile will not 
cause it to fail in tension, and the combined thermal deformation of the wing skin 
and the TPS tile is unlikely to cause tearing apart of the SIP layer. 


National Aeronautics and Space Administration 
Ames Research Center 
Drydcn Flight Research Facility 
Edwards, California, July 25, 19S6 
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APPENDIX — INTEGRATION OF EQUATION (10) 


Rewriting equation (10), 



The integration of this equation may be carried out as follows* Let 



Then, equation (40) becomes 



Equation (45) is the equation (11)* 



(48) 


Since h << p, the integrand of equation (40) may be expanded to give 

Pr •. /„u\2 


■ n ,+ *(?) a ~ 2 ? + •••]*■ 

Hsinq equation (41), equation (48) may be written as 

p -f l [’ *i(r ) 2 cos2 * + •"] d * 

- f [♦ - 1 (rf + 1 si " 2 ♦) * — ] 0 

-*[■♦* (f ) 2 ] 


or 


* - 1 + 




which is equation (13). 


i 


(49) 

(50) 

(51) 

(52) 
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